Axially symmetric static sources: A general framework and some analytical solutions 



CO 

o 

CN 

C 

c3 



L. HerreraQ A. Di Prisco and J. Ibanel] 

Departamento de Fisica Teorica e Historia de la Ciencia, Universidad del Pais Vasco, Bilbao, Spam 

J. Ospinci 

Departamento de Matemdtica Aplicada, Universidad de Salamanca, Salamanca, Spain. 

(Dated: January 14, 2013) 

We provide all basic equations and concepts required to carry out a general study on axially sym- 
metric static sources. The Einstein equations and the conservation equations are written down for 
a general anisotropic static fluid endowed with axial symmetry. The structure scalars are calculated 
and the inhomogeneity factors are identified. Finally some exact analytical solutions were found. 
One of these solutions describes an incompressible spheroid with isotropic pressure and becomes 
the well known interior Schwarzschild solution in the spherically symmetric limit, however it cannot 
be matched smoothly to any Weyl exterior metric. Another family of solutions was found that 
corresponds to an anisotropic fluid distribution and can in principle be matched to a Weyl exterior. 
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I. INTRODUCTION 

Observational evidence seems to suggest that de- 
viations from spherical symmetry in compact self- 
gravitating objects (white dwarfs, neutron stars), are 
likely to be incidental rather than basic features of these 
systems. This explains why spherical symmetry is so 
commonly assumed in the study of self-gravitating com- 
pact objects. 

However, the situation is not so simple. Indeed 
(putting aside the evident fact that astrophysical ob- 
jects are generally endowed with angular momentum, 
and therefore excluding all stationary sources), it is well 
known that the only regular static and asymptotically 
flat vacuum spacetime posesing a regular horizon is the 
Schwarzchild solution and all the others Weyl exte- 
rior solutions exhibit singularities in the curvature 
invariants (as the boundary of the source approaches the 
horizon). This in turn implies that, for very compact 
objects, a bifurcation appears between any finite pertur- 
bation of Schwarzschild spacetime and any Weyl solu- 
tion, even when the latter is characterized by parameters 
arbitrarily close to those corresponding to spherical sym- 
metry (see and references therein for a discussion 
on this point). 

From the above comments it should be clear that a 
rigorous description of static axially symmetric sources, 
including finding exact analytical solutions, is a praise- 
worthy endeavour. 

Accordingly, in this work we provide all ingredients 
and equations required for such a study. Of course, this 
issue has already been considered by several authors in 
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the past. Without pretending to be exhaustive in the 
revision of the literature on this problem, let us mention 
the pioneering paper by Hernandez Jr. [l2| . where a 
general method for obtaining solutions describing axially 
symmetric sources is presented. Such a method, or some 
of its modifications were used in [l3[ - [l5[ , to find sources 
of different Weyl spacetimes. 

This problem has also been considered in (l6j-[l^|. 
However in all these last references the line element has 
been assumed to satisfy the so called Weyl gauge, which 
of course severely restricts the family of possible sources 
(see the next section). 

In this work we present a general description of ax- 
ially symmetric sources by deploying all relevant equa- 
tions without resorting to the Weyl gauge, and consider- 
ing the most general matter content consistent with the 
symmetries of the problem. 

With this purpose in mind it would be useful to in- 
troduce the so called structure scalars. These form a set 
of scalar functions obtained from the orthogonal split- 
ting of the Riemann tensor. They were originally defined 
in the discussion about the structure and evolution of 
spherically symmetric fluid distributions. Such scalars 
(five in the spherically symmetric case) were shown to be 
endowed with distinct physical meaning [2p| - [23j . 

In particular they control inhomogeneities in the en- 
ergy density [20(, and the evolution of the expansion 
scalar and the shear tensor [2oM23| . Also in the static 
case all possible anisotropic solutions are determined by 
two structure scalars [201 ] . 

Furthermore, the role of electric charge and cosmolog- 
ical constant in structure scalars has also been recently 
investigated [24fj . 

More recently, such scalars and their applications were 
discussed also in the context of cylindrical [Hj], [26| and 
planar symmetry [27| . 

A set of differential equations for some of these scalars 
allows us to identify the inhomogeneity factors. 
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Finally we exhibit two families of solutions. One of 
them corresponds to an incompressible spheroid with 
isotropic pressure. It cannot be matched smoothly to any 
Weyl exterior spacetime. The second one corresponds to 
an anisotropic fluid and in principle is matchable to a 
Weyl exterior. 

Our paper is organized as follows: In the next sec- 
tion we shall describe the line element corresponding to 
the most general non-vacuum, axially symmetric static 
spacetime. Next we provide a full description of the 
source that is represented by a general anisotropic mat- 
ter. Einstein equations and conservation equations are 
explicitly written down for such a system. We also cal- 
culate the electric part of the Weyl tensor (its magnetic 
part vanishes) as well as the electric and magnetic part of 
the Riemann tensor. With this information we are able 
to obtain all the non-vanishing structure scalars corre- 
sponding to our problem. Two differential equations for 
such scalars allow us to identify the inhomogeneity fac- 
tors. The two families of solutions found are described 
in sections V and VI. A summary of the obtained results 
as well as a list of some unsolved issues are presented in 
section VII . Finally an appendix with the expressions for 
the components of the electric Weyl tensor, is included. 



II. THE METRIC AND THE SOURCE 

We shall consider bounded, static and axially symmet- 
ric sources. For such a system the most general line ele- 
ment may be written in cylindrical coordinates as: 

ds 2 = -A 2 dt 2 + B 2 [(dx 1 ) 2 + (dx 2 ) 2 } + D 2 dcj) 2 : (1) 

where A,B,D are positive functions of x 1 and x 2 . We 
number the coordinates x° — t, x 1 = p, x 2 = z, x 3 = tf>. 

We shall work in "Weyl spherical coordinates" (r, 9) 
defined by: 

p = rsm9, z = rcos6. (2) 

In these coordinates the line element reads: 

ds 2 = -A 2 dt 2 + B 2 (dr 2 + r 2 d0 2 ) + D 2 d<p 2 , (3) 

It is important not to confound these coordinates with 
Erez-Rosen coordinates (f, 9) given by 

p 2 = (f 2 — 2mf) sin 2 9, z = (r — m) cos 0, (4) 

where m is a constant to be identified with the monopole 
of the source. 

It should be stressed that our line element is defined 
by three independent functions, unlike the vacuum case 
where it is always possible to reduce the line element so 
that only two independent metric functions appear. In 
the interior this is not possible in general, though obvi- 
ously one may assume that as an additional restriction 



(the so called Weyl gauge), which amounts to assume 
that FLl + Rl = 0. 

In our notation the Weyl gauge is expressed by 
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r sin 9 
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Let us now provide a full description of the source. In 
order to give physical significance to the components of 
the energy momentum tensor, we shall apply the Bondi 
approach [28| . 

Thus, following Bondi, let us introduce purely locally 
Minkowski coordinates (r,x,y,z) (or equivalently, con- 
sider a tetrad field attached to such l.M.f.) by: 

dr = Adt; dx = Bdr; dy — Brd9; dz = Dd(j). 

(6) 

Denoting by a hat the components of the energy- 
momentum tensor in such locally defined coordinate sys- 
tem, we have that the matter content is given by 
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where p, P xy , P xx , P yy , P zz denote the energy density 
and different stresses, respectively, as measured by our 
locally defined Minkowskian observer. 

Also observe that P xy = P yx and, in general P xx ^ 

Pyy ^ Pzz- 
Introducing 

V a - (-1,0,0,0); K a = (0,1,0,0); L a = (0, 0, 1, 0), 

(8) 

we have 

T afi = (p + P zz )V a V p + P zz r) a0 + (P xx - P zz )K a Kfs 

+ (Pyy - Pzz)L a Lp + 2P xy K {a L fj) , (9) 

where r] a p denotes the Minkowski metric. 

Then transforming back to our coordinates, we obtain 
the components of the energy momentum tensor in terms 
of the physical variables as defined in the l.M.f. 

T a p = (p + P zz )V a V + P zz g af) + {P xx - P zz )K a K p 
+ [Pyy - Pzz)L a L l3 + 2P xy K {a L 0) , (10) 

where 

V a = {-A, 0, 0, 0); K a = (0, B, 0, 0); L a = (0, 0, Br, 0), 

(11) 

Alternatively we may write the energy momentum tensor 
in the "canonical" form: 



T a p = {p + P)V a V p + Pg afi + n a/S , 



(12) 



with 



U aP = (P xx - P zz ){K a K p - 

+ {Pyy ~ Pzz){L a Lp - ^) + 2P xv K (a Lp ) (13) 
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and 
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With the above information we can write the Einstein 
equations, which read: 
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where prime and subscript 9 denote derivatives with re- where the overdot denotes derivative with respect to t, 
spect to r and 9 respectively. and 

Also, the nonvanishing components of the conservation 
equations T a / = yield: 
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2P xy =0. (22) 



Equation (|20j) is a trivial consequence of the staticity, 
whereas (|2"Tj) and (l2"2"j) are the hydrostatic equilibrium 
equations. 



III. THE STRUCTURE SCALARS 



We shall calculate here the structure scalars for the 
static axially symmetric case. For that purpose, let us 
first obtain the electric part of the Weyl tensor (the mag- 
netic part vanishes identically). 
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The components of the electric Weyl tensor can be ited in the Appendix, 
obtained directly from its definition, 



Euu — CnavB V a V 13 , 



(23) 



Equivalently, the electric part of the Weyl tensor may 



where C^ av p denotes the Weyl tensor. These are exhib- also be written as: 
I 



1 



E a p = £\ (K a Lp + L a Kp) + £ 2 K a K/3 — -h a p + £ 3 [ L a Lp — -h a 



1 



(24) 



r 



where explicit expressions for the three scalars £\ , £2, £3 
are given in the Appendix. 

Next, let us calculate the electric part of the Riemann 
tensor (the magnetic part vanishes identically), which is 
defined by 



ip — V v *i a p^. 



After some lengthy calculations we find; 



(25) 



Y a p = Y TFl (K a L f:j + KpL a ) + Y TF2 ( K a Kp - -h a p 



Ytf 3 [ L a Lp — -h a p I + -Yrhap, 



(26) 



Finally, we shall find the tensor associated with the dou- 
ble dual of Riemann tensor, defined as: 



X, 
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K.psVV* = -r, a J P K p p S VW\ (31) 



with R* a a 5 = \r\ t(n (,R a n P ■ Thus, we find 



X a j3 — X T f 1 (K a Lp + KpL a ) + X T f 2 [ K a Kp — -h a p 



Xtf 3 [ L a Lp — -h a p 



-X T h a p, 



where 



Xt = 8717X, 
X TFl =-{£i+A-KP xy ), 



(32) 

(33) 
(34) 



where 



Y T = MP + P XX + Pyy + P„), (27) 



Y TFl =£ x - 4ttP xv , (28) 



Y TF2 =£2- MPxx - Pzz), (29) 



X T f 2 = -[£2+ MP XX - Pzz)] , (35) 



Xtf 3 = -[£3+ MPyy - Pzz )} ■ (36) 

The scalars It, Ytfi, Ytf2,Ytf3, Xt, Xtfi, Xtf2, 
Xtf3, are the structure scalars for our problem. 



IV. DIFFERENTIAL EQUATIONS FOR THE 
STRUCTURE SCALARS AND THE 
INHOMOGENEITY FACTORS 



Two differential equations for the Weyl tensor may 
be obtained using Bianchi identities pOj , they have 
been found before for the spherically symmetric and the 
cylindrically symmetric cases (see [25j , [31| and references 
Ytf 3 — £3 — 4>n(Pyy — Pzz)- (30) therein). Here we calculate them for our case. We obtain: 
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which, using (j2"T)) - ([3T))) and (|3"3"j) - l5()|) . may be written in 
terms of structure scalars: 



= _i [ Xtfw + X TF1 (ln B 2 D) ] - 



-X' TF2 +X TF2 (\nBD r y 
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Let us now turn to the inhomogeneity factors. The 
inhomogeneity factors (say *5>i) are the specific combina- 
tions of physical and geometric variables, such that their 
vanishing is a necessary and sufficient condition for the 
homogeneity of energy density (i.e. for the vanishing of 
all spatial derivatives of the energy density) . 

In the spherically symmetric case it has been shown 
that in the absence of dissipation the necessary and suf- 
ficient condition for the vanishing of the (invariantly de- 
fined) spatial derivative of the energy density is the van- 
ishing of the scalar associated with the trace free part of 

X a0 (see m, [H). 

We shall now identify the inhomogeneity factors in our 
case. 

First, observe that from (l39l) and (|40]) it follows at 
once that Xtfi = Xtf2 = Xtf3 = =>■ // = /ig = 0. In 
order to identify the above scalars as the inhomogeneity 
factors we need to prove that the inverse is also true (i.e. 

M = Me = =** Xtfi = X F F2 = Xtf3 = 0). 

For that purpose we shall first establish the behaviour 
of different variables in the neighborhood of r ~ 0. We 
shall demand that both A and B are regular functions, 
and 



r sin 6, 



(41) 



at r = 0. Then from (|2"Tj) it follows that in the neighbor- 
hood of r = 

and from (fT9l 



A e (0,9)^r [ \ B e (0,9)^r [ \ A(O,9)' ^r 2 



Using the above in (|A.9|) we have that 
£i(O,0)«r, 
implying because of (13"4")) 



(44) 



(45) 



Finally we shall assume that the three structure scalars 
Xtfi, Xtf2, Xtf3 are analytical functions (class C°°) in 
the neighborhood of r = 0. 

Then, if we assume \j! = [ig = 0. Evaluating (|4"0")) in 
the neighborhood of r ss 0, since X' TF1 (r = 0) is regular, 
it follows 



and 



X TF1 (r = 0) = 0, (46) 



XtF28 — 2XtF38 — XtF3 cot 9 = 0, (47) 



where (j4Tj) , ([43]) and (f46|) have been used (in what follows 
it is understood that all expressions are evaluated at r sa 
0). 

Next, from (l39l) we obtain 



*2Xt F2 — X 



TF3, 



(48) 



where (|4"6"| and the regularity of first derivatives of struc- 
ture scalars have been used. 

Then, feeding back (gSJ) into (g7|) we obtain 



TF2 



sin 



2/3, 



(49) 



(43) where a is a constant, implying Xtf2 = Xtf3 = 0. 
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Thus in the neigborhood of r pa we have Xtfi — 

XtF2 — XtF3 = 0. 
Next, taking r-derivative of (|40| . and evaluating at 
r pa we get Xj, F1 — 0. Continuing this process it 



(n) 
TF1 



for any n > 0. 



follows that X 

Also, from the r-derivative of ([39]) evaluated in the 
neighborhood of r pa it follows 



2X XTF2 — X FF3 . 



(50) 



Feeding this equation back into the r-derivative of (|4"0)) 
produces 



y 



/3 



sm 



2/3 , 



(51) 



where (3 is a constant, implying -X^ F2 = X^ F3 = 0. It is 
not difficult to see that this procedure can be continued to 
obtain, in the neighborhood of r pa 0, X, 



(n) _ v (n) _ 



in) 

Xj,p 3 — 0, for any n > 0. Therefore we can continue 
analytically their value at the center, from which we infer: 

Xtfi = Xtfi = Xtf3 = [/, = |itg = 0. (52) 

This last result allows us to identify the three structure 
scalars Xtfi, Xtf2, Xtf3 as the inhomogeneity factors. 
We shall next find some explicit analytical solutions. 
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V. THE INCOMPRESSIBLE, ISOTROPIC 
SPHEROID 

We shall now find an analytical solution correspond- 
ing to a bounded spheroid with isotropic pressures and 
homogeneous energy density. From the results of the pre- 
vious section, P5 ]) -([5D ]) and (JM])-®, it is evident that 
such a solution is also conformally flat. 

Thus let us assume P xx = P yy = P zz = P, P xy = 
and ji = jUo = constant 

For simplicity we shall assume the boundary surface S 
to be defined by the equation: 



r = n = constant. 



(53) 



Then to satisfy Darmois conditions (continuity of the first 
and second fundamental forms) we demand that all met- 
ric functions as well as r derivatives, to be continuous 
across S (see [I2])- Obviously 8 derivatives of A,B,D 
and A',B',D' are continuous too across E. 

From the above and IT6|) and (fT9|) it follows that 



P = 0. 



(54) 



Under the conditions above (1211) and (|22|) can be inte- 
grated to obtain: 



and 



A ' 



(55) 



(56) 



where £ is an arbitrary function of its argument. Using 
boundary conditions ([54]) in ([55)) (151)1) it follows that: 



a 



A(r\,6) — const. = — , 
Mo 



C, = constant. 



(57) 



Since, as mentioned before, our solution is conformally 
flat, then using £ i = £ 2 = £3 = and P xy = 0, in ([TT))) 
and (| A.6|) we obtain 
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Introducing the auxiliary function A(r, I 
A(r,6)=A(r,e)B(r,6), 



and assuming 



D( 



B(r,l 
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defined by 
(60) 

(61) 



the equations ([5"5)) and ([S^t can be integrated to obtain 



A(r,0)=A{r)+r X {0), 
R(r) + ruj{V) 



(62) 



where A, y, i? and lj are arbitrary functions of their 
argument. 

Next, from (jA.7j) and (|A.8|) . taking into account (f6T|) 
and ([6^]) we get 



X = a cos 6, A(r) = ar 2 + (3. 



(63) 



Where a, a and (3 are constants of integration. 
From the above it follows that conformally flat solu- 
tions are described by the line element 



ds 1 = 



1 



[R{r)+ruj{Q)Y 



[~(ar 2 + (3 + ar cos 6 fdt 2 + dr 2 + r 2 d6 2 + r 2 sin 2 (6»)d0 2 



(64) 



Next, from the condition P yy — P zz = and the equa- tions (fi~7 j) -([T8 j) we get 



cj(0) = bcos9, 



(65) 
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and from the condition P xx —Pyy 



and the equations (T16l) - (|T9)) we obtain 



R(r) = 7r 2 + 6, 



(66) 



where 6, 7 and 5 are constants of integration. 
Finally, the metric of incompressible conformally flat 
isotropic fluids can be written as follows. 



ds 2 = 



(■yr 2 + S + br cos 0) 2 



[-(ar 2 +f3 + ar cos 9) 2 dt 2 + dr 2 +r 2 d6 2 + 



r 2 sin 2 



(67) 



Next, the physical variables can be easily calculated. 
Thus, using ([57)1 into (|T5|) the energy density reads: 



(68) 



8irfi = 12jS- 3b 2 



To obtain the pressure we shall use (|55|) and ([57]). 
which produce 



8ttP = (36 2 - 12 7 <5) 
where 

C = Mo- 



f 



ar\ + f3 jr 2 + S + br cos 



ryrf + S ar 2 + j3 + ar cos ( 
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(69) 



(70) 



in order to satisfy the junction condition (|54|) . 

It may be instructive to recover the spherically sym- 
metric case (the interior Schwarzschild solution) . In this 
case we have a = b = 0. 

To see how this comes about, let us perform the trans- 
formation 

= 0, i=t, j> = (j), (71) 



6' 



where overbar denotes the usual Schwarzschild coordi- 
nates. Then, it is a simple matter to check that (|67| and 
(|69[) are identical to the well known expressions charac- 
terizing the interior Schwarzschild solution: 



3(1 
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^ _ 2m(r) j 1 /2 _ ^ _ 2Mjl/2 
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3(1 _ 2M)l/2 - (1 - 2w_(f)y/ 2 

where m(f), M and fx denote the mass function, the 
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total mass and the radius of the sphere respectively, and 
the following relationships are satisfied 
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and 



m(ri) = M, 
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At this point it is pertinent to ask the question: to 
what specific exterior spacetime can we match smoothly 
our solution? This is a relevant question, since there 
are as many different (physically distinguishable) Weyl 
solutions as there are different harmonic functions. 

The answer to the above question is the following: our 
solution cannot be matched to any Weyl exterior, even 
though it has a surface of vanishing pressure. This is 
so because the first fundamental form is not continuous 
across the boundary surface. 

Indeed from the continuity of gu and components 
at r = ri, we have: 



and 



A w (r u 6) 



A w (r i,9) 



jr 2 



: jrl 



(i + ari cos 
5 + bn cos ' 



5 + bri cos0, 



(78) 



(79) 



where Aw (f 1 , 0) denotes the ^Jgtt of any Weyl exterior 
solution (evaluated on the boundary surface). It is a 
simple matter to check that the two equations above can- 
not be satisfied unless a = b = which corresponds to 
the spherically symmetric case. This result is in agree- 
ment with theorems indicating that static, perfect fluid 
(isotropic in pressure) sources are spherical (see (32j | and 
references therein) 

Observe that the above result is a consequence of (|6T|) . 
Therefore in order to find matchable solutions we should 
relax this condition. In this later case, of course, nei- 
ther isotropy of pressure nor homogeneity of the energy 
density is preserved. 
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The remaining possibility is trying to match on a 
boundary surface given by the equation r — r\ (9) . How- 
ever this does not seem to solve the problem since in our 
case A(ri,9) — constant. 



VI. ANISOTROPIC INHOMOGENEOUS 
SPHEROIDS 

In order to find solutions that could be matched to a 
Weyl exterior, we shall relax the condition of isotropy of 
pressure and energy density homogeneity. 

The introduction of the pressure anisotropy is well mo- 
tivated from purely physical considerations. Indeed, it is 
well known that compact objects (white dwarfs and neu- 
tron stars) are endowed with strong magnetic fields (see 
[33|-[39[ and references therein). On the other hand it 
has been shown that the effect of such a magnetic field 
on a degenerate Fermi gas manifests itself through the 
appearance of a strong anisotropy due to the magnetic 
stresses (see [4(|-[43j and references therein). This in 

I 



turn may severely affect some important characteristics 
of compact objects (see [44|]-[46j and references therein). 
Besides the magnetic field, local anisotropy of pressure 
may be produced by a variety of physical phenomena 
(see [47j and references therein). 

Thus, let us assume P xy = £ \ = £3 = P yy — P zz = 0, 
although £2 ^ and P xx ^ Pyy Then from the equa- 
tionfUU) and (|A.9|) we obtain the equations §>8) and p^ . 

Next, introducing the auxiliary functions A(r,9) and 
R(r, 9) defined by 



A(r,6) = A(r,9)B(r,9)r, 
D(r,9) = R(r,9)B(r,9)r, 



the equation (f58| and ([59]) can be rewritten as: 



A 



R 'e 
R 



(80) 



(81) 



From Pyy - P zz = £ 3 = and the equations (|A.8|) (|A.11|) 

we get 



~R + ~R r ~ ~A> 



1 / Agg Ag Rg 



.4 



A R 



(82) 



R! A' B' 



Bg^R± 2 ~g_ _-_^ L _^ L) _ 
oi 7 + D + ~> + r ^A R B R +A B* ' ° ' ~ U ' 



Bi 



Ae 
A 



Bgg 
B 



(83) 



In order to find a simple solution that satifies the bound- 
ary condition, we choose: 



R(r, 9) = sin 9[b{r 1 )+j cos 6} 
then from the equations (|8"TT) - (|83I) we find 



A(r,8) = 



b(n) 



cos ( 



a(r), 



(84) 



(85) 



where a, 7, /? and a(r), b(r) are arbitrary constants and 
functions of integration respectively. a(r), a and f3 have 
the dimension of an inverse of length, whereas the others 
are dimensionless. 



B{r,0) 



(5r ?sin 9 — b(ri)cos9 + b{r) 



(86) 



For the above metric, Einstein equations yield the fol- 
lowing expressions for physical variables. 



8tt^ = {firT + b) 
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\ f3rT + b ~ {(3rT + b) 2 ~ {(3rT + b) 2 
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4(/3r + b') 4/3S 
f3rT + b + (3rT + b 



(87) 



9 



+ 



{l3rT + bf 
2a! 



f -2a' (pr + b') 3(/?r- 

\{/3rT + b)(aT + a) + JfirT 
6(^r + b') 4/3S 



h>\2 



3(/3A sin <9) 2 



-b) 2 {/3rT- 
2a(3A 2 sin 2 6 



bf 



1 



aT + a /3rT + b /3rT + b 
2a£ &i+47cos6>" 
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(89) 



where 

7 sin 2 

A = &i+7Cos#, r = — feicos^, b(r\) = bi 

S = Acos6»-7sin 2 6i. (90^ 

The equations above describe a wide class of solutions 
that can be matched to any specific Weyl metric by an 
appropriate choice of functions and constants of integra- 
tion. Furthermore physically reasonable models can be 
obtained at least from slight deviations from spherical 
symmetry. 

VII. CONCLUSIONS 

We have established the general framework to carry 
out a systematic analysis of general static axially sym- 
metric sources. By "general" we mean that the Weyl 
gauge was not assumed and the matter description is the 
most general compatible with axial symmetry and static- 
ity. 

Thus, we started from the most general line element 
and considered a general anisotropic fluid as source of the 
exterior Weyl spacetime. Relevant equations were then 
written down and structure scalars were calculated. 

We have seen that the three structure scalars associ- 
ated with the tracefree part of the tensor X a p define the 
inhomogeneity factors. 

We have found an exact analytical solution represent- 
ing a spheroid of isotropic pressure with homogeneous 
energy density. In the spherically symmetric limit our 
solution becomes the well known Schwarzschild interior 
solution. Such an interior cannot be matched (except in 
the spherically symmetric case) to any Weyl exterior. 

Matchable solutions can be found by relaxing the con- 
ditions of isotropy of pressure and density inhomogeneity. 
The anisotropy of pressure was also justified on physical 
grounds. An example was given in the last section. In 



order to study the physical relevance of nonsphericity in 
the structure of the source it is necessary to match the 
above metioned solution to a specific Weyl exterior (so 
that the arbitrary parameters of the source could be re- 
lated to the parameters of the exterior metric) , but such 
a task is beyond the scope of this paper. 

In general, solutions as the one presented here and 
many others found by either analytical or numerical pro- 
cedures could provide answers to important questions re- 
lated to stellar structure, namely: 

• What is the limit of compactness of a static axially 
symmetric source? 

• How is the above limit related to (influenced by) 
some exterior parameters such as the quadrupole 
moment of the source? 

• How are, intrinsically nonspherical physical vari- 
ables ( e.g. P xy ), related to multipole moments 
(higher than monopole)? 

• Are specific Weyl exteriors related to specific 
sources? 
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Appendix: Expression for the components of the 
electric Weyl tensor 



The nonvanishing components as calculated from (|23|) 
are: 
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These components are not independent since they satisfy For the three scalars £\, £2, £3 we obtain: 
the relationship: 
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Or, using Einstein equations we may also write: 
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